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An extension of the renormalized Jellium model which allows to study colloidal suspensions con- 
taining trivalent counterions is proposed. The theory is based on a modified Poisson-Boltzmann 
equation which incorporates the effects of counterion correlations near the colloidal surfaces using 
a new boundary condition. The renormalized charges, the counterion density profiles, and osmotic 
pressures can be easily calculated using the modified renormalized Jellium model. The results are 
compared with the ones obtained using the traditional Wigner-Seitz (WS) cell approximation also 
with a new boundary condition. We find that while the thermodynamic functions obtained within 
the renormalized Jellium model are in a good agreement with their WS counterpart, the effective 
charges predicted by the two theories can be significantly different. 

PACS numbers: 



I. INTRODUCTION 



Strongly acidic or basic groups are often used to stabi- 
lize colloidal suspensions against flocculation and precipi- 
tation. The resulting electrostatic repulsion prevents the 
like-charged particles from coming sufficiently near one 
another where a short-range van der Waals attraction can 
result in an irreversible association. It is well known that 
the Poisson-Boltzmann (PB) theory describes accurately 
the properties of such colloidal suspensions containing a 
symmetric 1:1 electrolyte For these systems the 

correlations between small, strongly hydrated ions are 
weak, making the mean- field PB equation quasi-exact. 
This is no longer true when the counterions are multi- 
valent. In this case, electrostatic correlations between 
the ions can no longer be neglected and the mean-field 
PB equation looses its validity. Although the PB equa- 
tion predicts very similar behaviors for suspensions con- 
taining cither monovalent or multivalent counterions jjj, 
both experiments and simulations show that this is not 
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correct. Many interesting phenomena inherent to col- 
loidal suspensions with multivalent counterions, such as 



colloidal charge reversal 



, like-charged colloidal at- 



traction jgl, or the reversal of the elcctrophorctic mobil- 
ity [jj] can not be captured by the simple mean field ap- 
proach. The description of these phenomena requires the 
use of more sophisticated and numerically more demand- 
ing approaches, such as the Monte Carlo simulations [lo| . 



the i nteg ral equations 
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, or the density functional thc- 



To account for the electrostatic correlations in suspen- 
sions containing multivalent counterions, while preserv- 
ing the simplicity of the PB theory, dos Santos et al. 
proposed a Wigner-Seitz (WS) model with a modified 
boundary condition at the colloidal surface QQ. The 



15| that the va- 



model is based on the Shklovskii's idea 
lidity of the mean-field theory can be extended if the cor- 
relations between the condensed counterions are properly 
taken into account. In the context of the PB equation, 
this can be done using a modified boundary condition 
which must be satisfied by the mean-field potential at 
the colloidal surface. Using this new boundary condi- 
tion, the PB equation can be solved inside the Wigncr- 
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Seitz cell to yield the ionic density profiles from which 
the effective colloidal charge can be calculated using the 
Alexander prescription 1 1 3L fl"3| - 

The WS cell model is based on the assumption that 
the liquid state structure of a colloidal suspension is not 
very different from a crystal one. This allows one to avoid 
the complicated many-body problem by considering only 
the electrostatic interaction between one colloidal par- 
ticle and its counterions and coions. As a further ap- 
proximation, the geometry of the WS cell is taken to 
match the colloidal one (e. g. spherical). The electro- 
neutrality condition is imposed by requiring that the elec- 
tric field vanishes at the cell boundary. Thus, within the 
cell model there no direct electrostatic interactions be- 
tween the different colloidal particles. Nevertheless, the 
effective colloidal charge obtained by matching the lin- 
ear and nonlinear solutions of the PB equation at the 
cell boundary is often used within the Derjaguin-Landau- 
Verwey-Overbeek (DLVO) potential [17[ to account for 
particle-particle interaction. 

The cell model should be particularly appropriate for 
suspensions with large volume fractions of colloidal par- 
ticles, while its validity for dilute systems is questionable. 
In this paper we propose a more liquid-state-like point of 
view, which is particularly appropriate for suspensions 
with low volume fraction of colloidal particles. 

Inside a suspension each colloidal particle feels the field 
produced by other particles and their condensed counteri- 
ons. If we neglect the correlations between the colloidal 
particles they, together with their condensed counteri- 
ons, can be thought to provide a uniform background in 
which free microions move. The effective colloidal charge 
and the background charge must be determined self- 
consistently, resulting in a Renormalized Jcllium Model 
(RJM) [L8|. This theory is particularly appropriate for 
calculating the effective colloidal charges because within 
the RJM, unlike in the cell model, the interaction po- 
tential between the colloidal particles has precisely the 
DLVO form The RJM has been successfully used 



to calculate the effective charges and the structural and 
thermodynamic properties of colloidal suspensions con- 
taining 1:1 electrolyte paper the RJM 
will be extended to account for the electrostatic correla- 
tions in suspensions with trivalent counterions. 

The paper is organized as follows. In section II, the 
general aspects of the model will be outlined. In section 
III we will describe the RJM with a modified boundary 
condition. In section IV the results of the theory will be 
presented. Finally, the conclusions and discussion will be 
given in section V. 

II. THE MODEL 

We consider a system of colloidal particles of radius a 
and (negative) charge —Zb are q (uniformly distributed on 
the surfaces), and counterions of radius r c and charge aq, 
(q is the charge of proton) inside an aqueous solution of 
volume V. We adopt a primitive model (PM) description 
in which the solvent is treated as a uniform continuum 
of dielectric constant e. The overall charge neutrality 
requires N c a — NZt, are = 0, where N c and N are the 
particle numbers of counterions and colloids, respectively 
The typical length scale that characterizes the system is 
the Bjerrum length, defined as A& = which is 7.2 A, 
in water at room temperature. In order to maintain the 
simplicity of the model, image charge effects which can 
become non-trivial for multivalent counterions 22|, are 



neglected within the primitive model approach adopted 
here — colloidal particle has the same dielectric constant 
as the solvent. 

Because of strong electrostatic interaction between the 
counterions and the colloidal particles, many of the coun- 
terions become condensed onto colloidal surface [l[ Q . A 
'complex' composed of one colloidal particles with a layer 
of its condensed counterions can then be regarded as a 
single entity carrying an effective charge Z e ff <C Zbare 
lj, |2|. The charge neutrality condition then becomes 
Pfa — Z e ffp = 0, where p = N/V, and pj is the number 
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density of free unassociatcd counterfoils. 

Traditionally colloidal suspensions have been modeled 
using a crystal-like approximation of a single macroion 
with its counterfoils inside a Wigner-Seitz cell, the ra- 
dius of which is determined by the volume fraction of 
colloidal particles. Although this picture is appropriate 
at large concentrations — when strong correlations be- 
tween the charged particles lead to a crystal-like ordering 
— it might not be appropriate for dilute suspension. To 
avoid the WS cell approximation we will instead use a 
liquid-state RJM. 

Suppose we fix one colloidal particle at the origin. 
Far from this particle the counterion density profile will 
have a Boltzmann-like form, Pf(r) = pfe~^ qa ^ r \ where 
4>{r) is the mean electrostatic potential. Note that pf 
refers only to free, uncondensed, counterions. The den- 
sity of other colloidal particles and of their condensed 
counterions provides a uniform neutralizing background 
Pback = pZback ■ The mean electrostatic potential satisfies 
the Jellium-Poisson-Boltzmann (JPB) equation: 
c 2 

a \ / a" 



VV(r) = — f 1 - e~ a ^ r A + ^^S(r - a), (1) 



where tf>(r) = j3qcj)(r) is the reduced potential, and 
k 2 = AirXbPfa 2 = AwXbpZbacka defines the inverse ef- 
fective Debye screening length. The self-consistency con- 
dition requires that the effective charge, calculated from 
the far field solution of this equation, and the background 
charge must have the same value, Zback = Zeff- We 
should note that the screening of the electrostatic po- 
tential in the far field is produced only by the free (un- 
condensed) ions. A similar behavior is implicit within in 
the cell model if one tries to define the effective colloidal 
charge, as is done within the Alexander prescription [r^ . 
However, this renormalization is less transparent within 
the cell model than within the Jcllium formalism. 

Although the above mean-field equation works very 
well for monovalent ions a = 1, it becomes a rather 
poor approximation when suspension contains multiva- 
lent counterions (a > 1). In such cases, the strong elec- 



trostatic correlations between the condensed counterions 
lead to significant deviations from the PB theory. The de- 
viations can be so strong that they qualitative modify the 
behavior of suspensions containing multivalent counteri- 
ons. In such suspensions, one finds that the counterion 
condensation can become so strong as to reverse the sign 
of the effective colloidal charge. Furthermore, addition 
of a multivalent electrolyte can result in attraction be- 
tween like-charged colloidal particles, thus destabilizing 
suspension against flocculation. 

To include the effects of the counterion correlations one 
can proceed in a number of different ways. One approach 
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is to use a weighted-density functional theory 
count for the corrections to the mean-field electrostatic 
potential. Another approach is to use the integral equa- 
tions theory 24|. All of these methods, however, have 



their own drawbacks and are significantly more compu- 
tationally demanding than the simple WS cell PB theory 
or the RJM. 

We note that the counterion correlations arc the 
strongest among the condensed counterions, since these 
ions are in the closest proximity of each other. Follow- 
ing Shklovskii we will, therefore, attempt to include the 
counterion correlations within the RJM using a modi- 
fied boundary condition at the colloidal surface. The 
condensed counterions will be treated as a strongly cor- 
related fluid — a concentrated quasi-two dimensional 
plasma. On the other hand, in the bulk the concentra- 
tion of counterions is quite small, so that the correlations 
can be neglected and the mean field approach is still suf- 
ficient. By matching the two regimes, the JPB equation 
is recovered, but with a new boundary condition at the 
colloidal surface. From now on, we will restrict our at- 
tention to trivalent counterions, a = 3. 

III. THE THEORY 

Strong electrostatic interactions between the colloidal 
particles and their counterions lead to counterion con- 
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dcnsation. The condensed ions are in thermodynamic 
equilibrium with the free ions of suspension. Close to 
the colloidal surface the counterion chemical potential is, 

f3p sc = ln(A 3 p sc ) + pn,. + /3qa<f>(a + r c ), (2) 

where A is the de Broglie thermal wavelength and p sc 
is the course-grained density of condensed counterions. 
The correlational chemical potential p c is given by that 
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of a two-dimensional one component plasma 

PUc = -1.65r + 2.6ir 1/4 - 0.26 log(r) - 1.95, (3) 



were T 



is the plasma parameter. Far 



2(a+r c ) 

away, in the bulk solution, the chemical potential is well 
approximated using the mean field electrostatic poten- 
tial, 



Pli b (r) = ln[A 3 Pf (r)} + (3qacj)(r). 



(4) 



The thermodynamic equilibrium between the condensed 
counterions and the free ions requires equality of Eqs. (2) 
and (4), from which follows 



p f (r) = p ace PPe e -P<rt<Kr)-+(°+r e )) m (5) 



This equation correctly describes the density profile after 
a short distance 8 from the colloidal surface. This cut-off 
distance delimits the region where the microion corre- 
lations arc important and the mean field approximation 
breaks down. However, since the range of counterion cor- 
relations is quite small and the PB density profile vary 
smoothly, Eq. (5) can be extrapolated all the way to 
the colloidal surface. This simplification then results in 
a new boundary condition for the JPB equation at the 
colloidal surface, 



Pf(a + r c ) = p sc e /3Mc 



(6) 



The value of p S c can be obtained using the strong cou 
pling theory 
find 



26j and coarse graining procedure 



131. We 



Z bare^b 



87r(3.701)(a + r c ) 4 



(7) 



The strong dependence of the coarse grained density on 
the colloidal radius is a direct consequence of the con- 
tact theorem which states that the difference between 
the contact and the bulk density of counterions is pro- 
portional to the square of the electric field [27j . Since the 



bulk counterion density is much lower than the counte- 
rion concentration at the colloidal surface, the theorem 
requires that the density near the colloidal surface scale 
as l/(o + r c ) 4 [13 



Together, Eqs. (3), (6) and (7) 
provide a new relation between the density profile at the 
colloidal surface and the bare colloidal charge Z^, are . This 
should be contrasted with the usual boundary condition 
for the PB equation, \ r=a = Zba ^ Xb which, however, 
does not capture the strong counterion condensation re- 
sulting from electrostatic correlations at short distance 
from the colloidal surface. 

The calculation of the electrostatic potential from Eq. 
(1) is now quite straight forward. Far from the colloidal 
particle the electrostatic potential has exactly the DLVO 
form, 



V>(r) 



Abe e 



(1 + na) r 



(8) 



Now suppose that we know Z e ff = Zb ac k 7 r c , a, and the 
volume fraction rj = 4ira 3 p/3. Eq. (8) then provides the 
electrostatic potential and the electric field far from the 
colloidal surface. Using these as the initial conditions, 
we can numerically integrate the JPB equation (1) to 
obtain the density p/(a + r c ). Numerical integration is 
performed using the usual Runge-Kutta algorithm. Eqs. 
(6) and (7) can then be used to calculate the bare col- 
loidal charge. In practice, of course, we know the bare 
charge and would like to calculate Z e f f . This can be eas- 
ily done by coupling the JPB solver with a root finding 
subroutine, such as the Newton-Raphson method. For 
each Z e ff there is a corresponding Zb are . The root find- 
ing subroutine allows us to efficiently search the values of 
Z e ff to find the one that corresponds to the given Z\, are . 
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FIG. 1: Effective charge as a function of the bare charge for 
suspensions with trivalent counterions. The volume fractions 
are a) T) = 0.125, b) r? = 0.01, c) r] = 0.005 and d) r? = 0.001. 
The solid lines have been calculated using the RJM, while 
the dashed lines have been obtained using the WS cell model 
. The difference between the two models diminishes as the 
volume fraction decreases. 

IV. RESULTS 

Fig. 1 compares the effective charge as a function of 
the bare charge for several colloidal concentrations, using 
the modified boundary condition within the RJM (solid 
curves) and within the WS cell model (dashed curves) 
13(. In all the calculations, the counterion and the col- 
loidal radii used are r c = 2 A and a = 100 A, respec- 
tively. 

Although both models predict similar qualitative be- 
haviors, there is a significant quantitative difference be- 
tween the two effective charges. The effective charges 
calculated using the RJM lie below the ones calculated 
in the WS cell approach. A similar result was found for 
the monovalent ions: apparently, the colloid- colloid cor- 
relations, implicit in the cell model, reduce the counte- 
rion condensation 28[. The difference, however, becomes 





FIG. 2: Density profiles calculated using the RJM (solid lines) 
and the cell model (dashed lines). The volume fractions are 
a) r? = 0.125, b) r? = 0.01 and c) r) = 0.001. The plasma 
parameter is F = 4.413. We see very small discrepancy, which 
is further reduced as the volume fraction is decreased. 

smaller as the colloidal concentration decreases, as is also 
observed in the case of 1:1 electrolyte. 

Unlike the monovalent situation, the effective charge in 
the case of trivalent counterions is not a monotonically 
increasing function of the bare charge, as can be clearly 
seen in Fig.l. Instead, after reaching a maximum, the ef- 
fective charge decreases with further increase of the bare 
charge. This general trend is in perfect agreement with 
the findings of the Monte Carlo simulations 13 1. 



In Fig. 2, the counterion density profiles calculated 
using both Jellium and WS cell models for a fixed bare 
charge Zb are \ b /a = 41.8, are displayed for different col- 
loidal concentrations. We conclude that the difference 
in the density profiles is even less pronounced when com- 
pared with the corresponding discrepancy in the effective 
charges calculated using the two models. It is important 
to remember that the density profiles obtained using the 
above theory are only valid after some distance from the 
colloidal surface. At short distances, the strong-coupling 
regime dominates over the mean-field. The present the- 
ory coarse-grains the whole near-field region into a mod- 
ified boundary condition for the JPB equation. 

Besides the effective charges and the ionic distribu- 
tions, the thermodynamic properties can also be easily 
calculated in the framework of the rcnormalizcd Jellium 
and the cell models. The osmotic pressure P within the 
RJM is a function of the bulk counterion concentration 
HQ: 



47rA b a 2 /3P = 



377A5 (naY 



(9) 



Employing Eq. (9), we have calculated the osmotic pres- 
sure as a function of the volume fraction for two fixed 



bare colloidal charges, Z bare =41.8 and Z b 



100.0, 



using both the WS cell model and the RJM. The results 
are shown in Fig. 3. Again, we see only a very small dif- 
ference at moderate volume fractions. For small volume 
fractions, the osmotic pressures are identical within the 
two models. 

V. SUMMARY AND CONCLUSIONS 

We have extended the range of applicability of the 
rcnormalizcd Jellium model to describe suspensions with 
multivalent counterions in the absence of added salt. The 
model uses Shklovskii's idea to include the counterion 
correlations as a modified boundary condition for the 
JPB equation. A similar strategy has already proven 
to be successful for the modified WS PB model [3, . 




FIG. 3: Osmotic pressure as a function of the volume fraction 
for colloidal suspensions with trivalent counterions and bare 
charges Zbare = 41.8 (upper curves) and Zbare = 100.0 (lower 
curves). The corresponding plasma parameters are F = 4.413 
and r = 6.826, respectively. The solid lines have been calcu- 
lated using the RJM, and the dashed lines have been obtained 
using the WS cell model. 

Comparing the predictions of the renormalizcd Jellium 
model and the WS model, we find a quantitative dif- 
ference in the values of the effective charge of colloidal 
particles. Since the far field potential within the renor- 
malized Jellium formalism has precisely the DLVO form, 
we expect that the effective charges calculated using this 
formalism should be more reliable for structural calcula- 
tions. On the other hand the measurable thermodynamic 
quantities such as the osmotic pressure come out to be 
practically identical in the two models. 

Besides providing an alternative to the cell model, 
the RJM can be further extended to take into account 
the colloid-colloid correlations. The homogeneous back- 
ground charge distribution can be replaced by a non- 
homogeneous one related to the colloid distribution func- 
tion. This approach has been successfully implemented 
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Un- 



for suspensions with monovalent electrolyte 
fortunately, for the case of multivalent counterions, inclu- 
sion of colloidal correlations is not so straightforward. To 
have an accurate structure function one needs to know 
the interaction potential not only in the far field, but also 
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in the near field. The DLVO effective interaction poten- 
tial, however, is not valid at short distances where it be- 
comes^ strongly modified by the counterion correlations 
The work in this direction is now in progress. 
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